Fixed points of Higher Derivative Gravity 



O 

o 

(N 



(N 
> 
00 

o 
o 



Oh 

43 ■ 



X 



Alessandro CodelloEl 

Dipartimento di Fisica Teorica, Universita di Trieste, Viale Miramare, 1-34014 Trieste, Italy 

Roberto Percaccfj] 

SISSA, via Beirut 4, 1-34014 Trieste, Italy, and INFN, Sezione di Trieste, Italy 

We recalculate the beta functions of higher derivative gravity in four dimensions using the one- 
loop approximation to an Exact Renormalization Group Equation. We reproduce the beta functions 
of the dimensionless couplings that were known in the literature but we find new terms for the beta 
functions of Newton's constant and of the cosmological constant. As a result, the theory appears to 
be asymptotically safe at a non-Gaussian Fixed Point, rather than perturbatively renormalizable 
and asymptotically free. 
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The earliest attempts at constructing a Quantum Field 
Theory (QFT) of gravity were based on the application 
of perturbative methods to Einstein's theory. It was soon 
understood that such methods would not succeed due to 
the perturbative nonrenormalizability of Einstein's the- 
ory. It was then natural to try with more general types 
of dynamics. Lagrangians with four derivatives of the 
fields give propagators that fall off with the fourth power 
of momentum, leading to improved convergence of loop 
integrals. It was indeed proven that a generalization of 
Einstein's theory containing terms quadratic in the cur- 
vature tensor is renormalizable in flat space perturbation 
theory [1]. It was also established in a series of papers 
[2-4] that the dimensionless couplings of this theory (the 
inverse coefficients of the curvature squared terms) are 
asymptotically free. The beta functions of the dimen- 
sionful couplings - Newton's constant G and the cosmo- 
logical constant A - are gauge-dependent, but the beta 
function of the dimensionless product AG is not, and this 
variable has also been claimed to be asymptotically free, 
justifying the use of flat space perturbation theory. 

A perturbatively renormalizable and asymptotically 
free QFT holds to arbitrarily high energy scales, so this 
could be regarded as a serious candidate for a fundamen- 
tal theory of quantum gravity. Unfortunately, it is not 
free of problems, the most notorious one being the ap- 
parent lack of unitarity: the "bare" action contains mas- 
sive negative-norm states (ghosts) at tree level. It was 
pointed out in [2,6] that these ghosts may not correspond 
to physical particles when quantum effects are taken into 
account, but there exists to date no convincing proof that 
this happens. Another, less well-known problem is that 
asymptotic freedom of AG requires the choice of the un- 
stable fixed point -5.467 for the parameter uo of eq.(l) 
below, see [5]. Pending progress on these issues, higher 
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derivative gravity does not seem to have gained wide ac- 
ceptance as a fundamental theory. For a review of higher 
derivative gravity see [5]; for the state of the art see [7]. 

Nowadays, it is understood that Einstein's theory and 
its higher derivative generalizations can be successfully 
treated as effective QFTs with a cutoff presumably close 
to the Planck scale [8]. This is more than enough to 
cover all available experimental data, but this fact has 
not stopped the search for an "ultraviolet completion" of 
the theory. From the preceding remarks, it would seem 
that QFT can only work if some nonperturbative mech- 
anism is invoked, the most promising one being as fol- 
lows. A QFT that admits a Fixed Point (FP) with a 
finite number of UV-attractive directions can be predic- 
tive and hold to arbitrarily high energies. This behaviour 
was called "asymptotic safety" in [9]. A perturbatively 
renormalizable and asymptotically free theory is a spe- 
cial case of asymptotically safe theory, where the FP is 
the Gaussian FP (a free theory). More general asymp- 
totically safe theories will be based on nontrivial FPs. 

The question arises whether a QFT of gravity could 
have this behavior. The first positive evidence came long 
ago from studies in 2 + e dimensions [9,10], but techni- 
cal issues then slowed down progress on this front for 
some time. In the last ten years, using an Exact Renor- 
malization Group Equation (ERGE), the existence of a 
nontrivial FP has been established in four dimensions 
for a truncation of the action containing the cosmologi- 
cal and Einstein-Hilbert terms [11,12], also in the pres- 
ence of matter fields [13,14]. Independent evidence for a 
nontrivial FP also comes from Monte Carlo simulations 
[15,16]. However, so far only partial results are known 
for higher-derivative terms [17,18,14]. 

The behaviour of A and G in this approach is quite dif- 
ferent from the one predicted in the literature on higher 
derivative gravity. In order to make a direct comparison, 
we have recalculated the beta functions of higher deriva- 
tive gravity, starting from a one-loop approximation of 
the ERGE. We find some important modifications in the 
beta functions of Newton's constant and of the cosmolog- 
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ical constant, in such a way that the theory appears to 
be asymptotically safe at a nontrivial FP, rather than at 
the Gaussian FP. We report here the main results; details 
will be given elsewhere. 

A general (Euclidean) theory containing terms 
quadratic in curvature has an action of the form 
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where Z = l/l6irG, C 2 is the square of Weyl's tensor, 
E is the integrand in Euler's topological invariant x = 
/ dx^fgE. We neglect the total derivative V 2 i?. 

For a quantum treatment, this action has to be sup- 
plemented by the gauge-fixing term, which is chosen to 
be of the form 



Sgf = J Jxy/gxpY^Xt 



(2) 



where Xv — V^/i M „ + /3 V ' v h (all covariant derivatives 
are with respect to the background metric) and Y^ v = 
i [g^ V 2 + 7VV" - . The ghost action con- 

tains the term 



Sr. 
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where (A gh )» 
a term 



-5 V U - (1 + 2/3) V M V + R v as well as 



S h = - 



l - Jd 4 x^b^b p 
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due to the fact that the gauge averaging operator Y de- 
pends nontrivially on the metric. We follow earlier au- 
thors in choosing the gauge fixing parameters a, /3, 7 and 
S in such a way that the quadratic part of the action is: 



(r fc + S GF ) 



(2) 
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where A' 4 ' = ID 2 + V A V P V A + U. For details of the 
operators K, V and U we refer the reader to [7], whose 
notation we mostly follow. 

The main tool in deriving nonperturbative information 
about the theory is the gravitational ERGE [19] 
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where is a coarse-grained effective action depending 
on a momentum scale k and the kernels Rfc act as infrared 
cutoffs. 

In order to derive the beta functions of the couplings 
A = k~ 2 A, G = k 2 G, A, uj and 6, we assume for Tfc the 
form (1) and insert it, together with the gauge-fixing and 
ghost terms (2,3,4), into the ERGE. Then, to calculate 
the r.h.s. of the ERGE we choose the cutoffs as follows: 



lR^(A (gh) ), 

(Y) = 1R (2) {Y), where R^\z) is a suitable pro- 
file function chosen to suppress the propagation of field 
modes with momenta below k. We will use the so-called 
optimized cutoff [20] i? fc n) '(z) = (ak n - z)9{ak n - z), with 
a = 1 unless otherwise stated. 

We restrict ourselves to the one-loop approximation, 
which in the context of the ERGE consists of taking into 
account only the explicit dependence of Rk(z) on k, ne- 
glecting the implicit dependence due to the presence of 
running couplings in the cutoff function. (In the case 
of the Einstein-Hilbert action, where the r.h.s. of the 
ERGE can be computed exactly, it is known that this 
approximation does not change the general behaviour.) 



The traces are evaluated with heat kernel methods, keep- 
ing all terms up to B4, and using the results of [21]. This 
procedure provides a logically and computationally inde- 
pendent derivation of the beta functions. 

The beta functions of the dimensionless couplings ap- 
pearing in (1) turn out to be: 
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They agree with those calculated in dimensional regular- 
ization [4,5,7]. The coupling A has the usual logarithmic 
approach to asymptotic freedom, while the other two cou- 
plings have the FP values w (-5.467, -0.0228) and 
6* w 0.327. Of the two roots for u, the first turns out to 
be UV-repulsive, so the second has to be chosen [4,5,7]. 
The beta functions of A and G are: 
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PIG. 1: The flow in the (A, G)-plane 



where q(u) = (83 + 70a; + 8o; 2 )/187r. The first two terms 
in each beta function exactly reproduce the results of 
[4,5,7], the remaining ones are new. The origin of the 
new terms can be easily understood. The beta func- 
tions were originally derived as coefficients of 1/e poles 
in dimensional regularization, which correspond to loga- 
rithmic divergences in the effective action. In our heat 
kernel derivation these terms are given by the B4 coeffi- 
cient. The new terms that we find come from the B2 and 
-Bo coefficients, which in a conventional calculation of the 
effective action would correspond to quadratic and quar- 
tic divergences. Dimensional regularization is ill-suited 
to compute these terms. It is important to stress that 
our "Wilsonian" calculation of the beta functions does 
not require any UV regularization. The only ambiguity 
is in the choice of the cutoff functions, but no reasonable 
choice could remove the B 2 and B terms. 

To picture the flow of A and G, we set the remaining 



variables to their FP values uj = a;* 



= 8 t , and A 



A* = 0. Then, defining = q{u)*) « 1.440 the flow 
equations (8) can be solved analytically: 



A(t) 



G(t) 



(27rA -G (l- e 4t ))e 
7r(2- g *G (l-e 2 *)) 
2G e 2t 
2 - q,G (l - e 2t ) ' 



-2/ 



(9a) 
(9b) 



The resulting flow in the (A, G)-plane is shown in 
Fig.l. It has two FPs: the Gaussian FP at A = G = 



and another one at 
1 
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The attractivity properties of these FPs are determined 
by the stability matrix 



Mi. 



-2 - q*G I - g*A 
2 - 2q t G 



At the Gaussian FP the eigenvalues of M are (—2, 2); the 
attractive eigenvector points along the A axis and the re- 
pulsive eigenvector has components (1,2tt). As expected 
on general grounds [9], the eigenvalues are the opposite 
of the canonical dimensions of A and G. At the non- 
Gaussian FP the eigenvalues of M are (—4, —2) with the 
same eigenvectors as before. The FP given by (10) is UV- 
attractive in all five couplings. Note the "critical" trajec- 
tory joining the Gaussian to the non-Gaussian FP, which 
is tangent to the repulsive eigenvector in the origin and 
is actually given by G(t) = 27rA(t) for all — 00 < t < 00. 

From this calculation one can derive some physical pre- 
dictions. The first is the UV-limit of the cosmological 
constant in Planck units AG = AG, which is well known 
to be gauge-independent and is also independent of the 
cutoff parameter a. In contrast to [4,5,7], we find that 
AG tends to the finite value 2/(7rq 2 ) f=a 0.307. Of course 
this is an asymptotic UV value and to compare it with 
cosmological observations one would have to run the RG 
down to extremely low values of k. 

Another prediction is the asymptotic value — 2o;*/3 w 
0.0152 for the ratio between the coefficients of B 2 and 
C 2 . It is interesting to observe that the flow induced 
by a large number N of minimally coupled matter fields 
gives for this ratio the value 5ns/(3ns + I8nn + 3Qum), 
where 715, nr> and nu are the numbers of scalar, Dirac 
and gauge fields [14]. This number is also quite small in 
realistic unified theories. Thus both with and without 
matter it seems that, in the UV limit, fluctuations of the 
conformal tensor will be more suppressed than fluctua- 
tions of the Ricci tensor. 

The flow that we find here is almost identical to the 
flow obtained in the large N limit [14], where the coeffi- 
cients u>*, 9„ and g* are determined by ns, no and um- 
A remarkable feature of the large TV limit, in conjunc- 
tion with the use of optimized cutoffs, is that all higher 
powers of curvature are absent at the FP. This raises the 
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hope that asymptotically safe gravity may be describable 
by a finite number of terms in the action (generically, 
one would expect to have infinitely many terms, with re- 
lations between the coefficients such that only a finite 
number of parameters is left arbitrary) . 

Our flow is also similar to the one obtained in the 
Einstein-Hilbert truncation [12], where, however, the 
critical exponents at the non-Gaussian FP are complex, 
resulting in a spiralling approach to the FP. This similar- 
ity may be somewhat surprising, because in the Einstein- 
Hilbert truncation the higher derivative terms are absent 
while here they dominate the dynamics. To some extent 
it can be understood by the following argument. In grav- 
ity at low energies the couplings do not run, and therefore 
the relative importance of the terms in the action can be 
determined simply by counting the number of derivatives 
of the metric. For example, at low momenta p <C \[Z (re- 
call that Z is the square of the Planck mass) , the terms 
in the action (1) with four derivatives are suppressed rel- 
ative to the term with two derivatives by a factor p 2 jZ. 
This is not the case in the FP regime: if we consider 
phenomena occurring at an energy scale p, then also the 
couplings should be evaluated at k « p. If there is a 
nontrivial FP, Z runs exactly as p 2 and therefore both 
terms are of order p A . This is just a restatement of the 
fact that in the Einstein-Hilbert truncation the graviton 
has an anomalous dimension equal to two, making its 
propagator behave like p~ A at high energy. 

Partial results for the four-derivative couplings, but 
going beyond one loop, have been derived using the 
ERGE in [17]. Using a spherical background, where 
/ d 4 x^/gC 2 = 0, / d 4 x^/gR 2 = 384tt 2 and x = 2, the 
beta function of the combination — ^ + f can be 
derived. In the absence of further input it is impossible 
to disentangle the beta functions of the individual cou- 
plings. Nevertheless, this provides valuable information. 
In particular, since a finite FP-value is found for a com- 
bination of couplings, this calculation suggests that the 
asymptotic freedom of A, X/ui and X/8 that we find here 
may be only a consequence of the approximations that 
we made, and that in a more accurate calculation some 
or all of these coefficients will reach finite values instead 
of running logarithmically. One also expects, as in [17], 
that the degeneracy of the stability matrix is lifted and 
that all couplings are either relevant or irrelevant. 

To summarize, we have found that within our approx- 
imations higher derivative gravity has a fixed point with 
the following properties: A and G are nonzero and UV- 
relevant, while the couplings of the terms quadratic in 
curvature are asymptotically free and marginal. Experi- 
ence with the Einstein-Hilbert truncation suggests that 
the FP will persist in a more precise treatment, up to a 
finite shift of the FP-values of the couplings, and of the 
critical exponents. The Gaussian FP is unstable: even 
an infinitesimal value for G will generate a nonvanishing 
A and push the system towards the other FP. 



Among other things, these results solve the second of 
the problems mentioned in the introduction. Concerning 
the issue of unitarity, we can say, from our Wilsonian 
point of view, that the presence of ghost poles at the 
Planck scale has to be assessed by considering the ac- 
tion Tfe for k f=a mpi anc k, which is probably quite different 
from the FP action. Thus, tree level analyses of the FP 
action are of little significance, as already pointed out in 
[2,3,6]. This is generally accepted in the case of QCD: 
a tree level analysis of the QCD FP action would pre- 
dict the existence of states that are not observed in the 
physical spectrum, but this is no longer considered a seri- 
ous argument against this theory. In view of this, and of 
the results reported here, we think that higher derivative 
gravity deserves renewed attention. 
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